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We investigate medium effects in QCD like chromoelectric screening and quasi-particle 
mass generation by calculating the heavy quark potential as well as the temporal quark 
and gluon Coulomb gauge propagators in quenched approximation. 

m ■ 1. INTRODUCTION 

O ■ 

i The perturbative description of the Quark-Gluon Plasma (QGP) breaks down even at 
very high temperatures |T]]. The reason for this is the generation of different length (time) 
scales at finite temperature. At high temperatures (T ^ ^qcd) the gauge coupling con- 
stant is small g <^ I and well separated length scales exist: 1/T <^ l/{gT) <^ l/{g'^T). 
For such high temperatures the region of applicability of the perturbation theory is es- 
Q^l tablished. While for distances smaller than 1/T the ordinary perturbation theory should 
be applicable, at distances smaller than 1/gT the hard thermal loop resummed perturba- 
tion theory should be applied. For distances larger than 1/g'^T, however, perturbation 
- theory is no longer valid [0. For interesting temperatures up to a few times Tc, however, 
', g>l and it is not clear up to which distances (times) the perturbation theory is appli- 
ed I cable. At short distances r -C one expects that medium effects are not important. 
For such small distances one should also expect the ordinary perturbation theory to be 
viable, because at zero temperature perturbation theory works at short distances. Nat- 
urally the question arises: at what distances the medium effects become important and 
up to which distance perturbation theory is applicable? In the present work we try to 
answer these questions by studying the heavy quark potential and temporal quark and 
gluon propagators in Coulomb gauge. 

2. THE HEAVY QUARK POTENTIAL 

We consider the color averaged heavy quark potential defined by 

■< L{R)L^{0) >- 



V{R,T) = -Tin 



<L> 



(1) 



Here L{R) = TtH^^^^ Uo{R,xo) is the Polyakov loop (A^^ is the temporal extent of the 
lattice). So far most studies have concentrated on the long distance behavior of the heavy 
quark potential [RT > 1). For the physics of heavy quarkonia at finite temperature, 
however, it is important to know the potential for RT < 1. In order to explore the short 
distance behavior of the heavy quark potential we performed simulations on 64^ x 16 
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lattices with the standard Wilson action at two different temperatures T = 1.5Tc and 
T = 3Tc, as well as on 32^ x 8 lattices with a tree level Symanzik improved action for 
r/Tc = 1.05, 1.14, 1.23, 1.50, 3.0, 6.0, 12.0. The color averaged potentialQ is related to 
the singlet (Vi)and octet (Vg) heavy quark potentials by the well-known formula 

V(R.T) = -Tl„(lexp(-15i^)) + ?exp(-l^)). (2) 

At very short distances RT <C 1 one would expect that the singlet and octet potentials 
are given by their perturbative zero temperature expressions Vi = Cia/R, with ci = 
—4/3 for the singlet and Cg = +1/6 for the octet case. Thus for small distances one 
would expect V{R,T) ~ 1/_R. Expanding Eq. (||) up to one recovers the well known 
perturbative result for the color averag ed potential V{R,T)/T = -aV(3i?T)' |]. We 
note, however, that the expansion of Eq. (H) couples the perturbative expansion to the 
high temperature expansion and therefore the above formula is obviously not valid at very 
short distances where |Vi/T| ^ 1. Motivated by the leading order perturbative result we 
define an effective couphng alff{R,T) = -9V{R,T)T[1/ R]~^. Here [1/R] denotes the 
lattice Coulomb potential. Our results on aeff{R,T) are summarized in Fig. la. One 
can see that there is a narrow region in RT where is approximately constant as 
one would expect from the leading order perturbative result. We also note that in this 
region of RT aeff{R,T) decreases with temperature as one would expect for the running 
coupling constant. For temperatures T < l.ST^ the effective coupling seems to decrease at 
small distances RT<0.2 which may imply that in this region the potential behaves like 1/R 
according to the arguments given above. For higher temperatures the 1/i? behavior sets in 
at smaller distances because the coupling a is smaller and the high temperature expansion 
works well even at shorter distances. This probably is the reason that this behavior 
so far is not seen in our data for T > 3Tc which are restricted by the lattice spacing 
a = 0.0625r-i for A^^ = 16 and 0.125r-i for = 8. For distances RT > 0.3 the effective 
coupling decreases with increasing RT which implies the onset of screening. Thus we can 
distinguish three different regions of RT which are characterized by qualitatively different 
behavior of V{R,T): the "true" short distance region where V ~ 1/-R, the intermediate 
perturbative region V = —a'^/{3RT)'^ and the "large" distance region where screening 
sets in. We note that for the lowest temperature T = l.OSTc the behavior of V{R,T) is 
compatible with 1/i?^ for RT > 0.3 which may imply that the effect of screening is very 
small. To quantify the effects of screening at higher temperatures it is useful to define 
an effective screening mass /ie//(-R, T) = —\n{—9V{R,T)[l/R]~'^/a'^(T))/R with a{T) = 
aeff{R = 0.125T~^,T). The behavior of fieff{R,T) as function of R is shown in Fig. 
lb. As one can see from the figure the effective screening masses increase with increasing 
distance but up to distances RT = 1.5 do not reach a plateau which can be identified with 
the true screening mass. In fact this behavior of the effective screening masses could be 
explained qualitatively in perturbation theory if one takes into account the momentum 
dependence of the gluon self energy Iloo{k) [^]. Note that the local screening masses in Fig. 
lb are always smaller than screening masses extracted from the large distance behavior 
of the heavy quark potential in El yU ~ 3T. They are also smaller than the perturbative 



^ We stick here to the commonly used notion of a potential, although for T > one actually calculates 
the heavy quark free energy. 
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Figure 1. The effective coupling and tlie effective screening masses at different tempera- 
tures. Tfie open symbols in tlie left figure are tlie results from simulations on a 64^ x 16 
lattice with the Wilson action, others correspond to the data obtained on a 32'^ x 8 lattice 
with an improved action. 

Debye screening mass g{T)T if g is defined by g{T) = ^A7caeff{R = 0.125T~i, T). We 
thus conclude that the color averaged heavy quark potential has a complex structure for 
interesting distances, RT < 1, and temperatures a few times Tc. 



3. THE TEMPORAL QUARK AND GLUON PROPAGATORS 

The temporal quark and gluon propagators are related to properties of quasi-particles 
in the QGP. In fact, in hard thermal loop (HTL) approximation these propagators are 
dominated by quasi-particle poles. [0. The quasi-particle picture of QGP finds appli- 
cation in resummed perturbative calculations of different thermodynamic quantities [0. 
We have performed simulations in quenched QCD with Wilson fermions on a 64^ x 16 
lattice at T = ST^.. We have used the standard Wilson action for the gauge fields and an 
0{a) improved fermion action (clover action). We have calculated the quark and gluon 
propagator Di{T,p) {i = q,g) in mixed (t, p) representation ( here p is the absolute value 
of the spatial momentum ). 

The propagator in this mixed representation can be related to the spectral function by 
= du p{lj , p) exp{—u!r) / {l±exp{—LJT)) . Here (i = q,g) and +(— ) correspond 
to the case of quarks (gluons). From the Monte-Carlo data for Di{T,p) we can extract 
the spectral functions p{uj,p) with the help of the Maximal Entropy Method (MEM) [Q. 
Since the quark and gluon propagator are gauge dependent quantities one has to fix a 
gauge. We have fixed the Coulomb gauge because in this gauge p{uj,p) is positive which is 
necessary for MEM. The peaks of the spectral function p{uj,p) define the dispersion rela- 
tion uj{p) which is expected to be gauge invariant. Since the smallest non-zero momentum 
Pmin = 1.57T is relatively large no detailed information about the properties of collective 
excitations like the longitudinal excitation in the gluon sector or the plasmino excitation 
in quark sector can be provided (these excitations exist only in the small momentum 
region [^). Here we will not extract the spectral function but postpone that for future 
publications. Instead we compare our data for temporal quark and gluon propagators with 
the predictions of perturbation theory in HTL approximation. In Fig. 2a we show our 
results for the transverse gluon propagator Dj(r, p) at momenta calculated on 40 gauge 
fixed configurations. We also show there the prediction of the HTL approximation using 
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Figure 2. The transverse gluon propagator (a) and the quark propagators (b) for 
different momenta. The sohd, dashed and dashed-dotted hues correspond to the prediction 
of perturbation theory in HTL approximation for momenta p/T = 0, 1.57 and 3.14 
respectively. 

a couphng constant ~ 1.6 suggested by the short distance behavior of the heavy-quark 
potential. As one can see from the figure the data disagree substantially with the HTL 
prediction even at momenta p = 3.14T. We note that corrections to the HTL approxima- 
tion will not resolve this discrepancy since they lead to smaller a;(p) values shifting the 
propagator to values larger than the HTL result 0. 

The quark propagator in mixed representation can be written as Dq{T,p) = 7oF(r, p) + 
^ ■ nG{T, p) + H{t, p) with n = p/p. In the chiral limit (m, = 0) H{t,p) = 0. For the value 
of the hopping parameter k, = 0.1339 (which corresponds to the quark mass close to the 
chiral limit) used in our simulation we have found H{t,p) to be zero within the statistical 
accuracy reached in our simulations. In Fig. 2b we show our results for F{t,p) calculated 
on 40 configurations and compared with predictions of the HTL approximation. As in 
the case of the gluon propagator we find large deviations from the HTL predictions. The 
situation is similar for G{t,p). We have found that there is no choice of g which can 
provide agreement between lattice results and HTL. We thus conclude that at T = 3Tc 
medium effects in temporal quark and gluon propagators are stronger than suggested by 
HTL perturbation theory. 
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